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I Abstract. The observation of neutrinos from Supernova 1987A has con- 

firmed the theoretical conjecture that these particles play a crucial role 
during the collapse of the core of a massive star. Only one per cent of 
the energy they carry away from the newly formed neutron star may ac- 
count for all the kinetic and electromagnetic energy responsible for the 
spectacular display of the supernova explosion. However, the neutrinos 
emitted from the collapsed stellar core at the center of the explosion cou- 
^ ■ pie so weakly to the surrounding matter that convective processes behind 

CN i the supernova shock and/or inside the nascent neutron star might be 

required to increase the efficiency of the energy transfer to the stellar 
mantle and envelope. The conditions for a successful explosion by the 
■ neutrino-heating mechanism and the possible importance of convection 

, in and around the neutron star are shortly discussed. Neutrino-driven 

Q^ ' explosions turn out to be very sensitive to the parameters describing 

^ , the neutrino emission of the proto-neutron star and to the details of the 

Qh' dynamical processes in the collapsed stellar core. Therefore uniform ex- 

Q . plosions with a well defined energy seem unlikely and type-II supernova 

^ I explosions do not offer promising perspectives for being useful as standard 

c/3 . candles. 
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Even after ten exciting years SN 1987A keeps rapidly evolving and develops new, 
unexpected sides like an aging character. In the first few months the historical 
detection of 24 neutrinos in the underground facilities of the Kamiokande, 1MB, 
and Baksan laboratories caused hectic activity among scientists from very dif- 
ferent fields. In the subsequent years the scene was dominated by the rise and 
slow decay of light emission in all wavelengths which followed the outbreak of 
the supernova shock and contained a fiood of data about the structure of the 
progenitor star and the dynamics of the explosion. Now that the direct emission 
has settled down to a rather low level, the supernova light which is reflected from 
circumstellar structures provides insight into the progenitor's evolution. Even 
more information about the latter can be expected when the supernova shock 
hits the inner ring in a few years. 

The neutrino detections in connection with SN 1987A were the final proof 
that neutrinos take up the bulk of the energy during stellar core collapse and 
neutron star formation. Lightcurve and spectra of SN 1987A bear clear evidence 
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Figure 1. Sketch of the post-collapse stellar core during the neutrino 

heating and shock revival phase. At the center, the neutrino emitting 
proto-neutron star (PNS) with radius R^s is shown. The average ra- 
dius where neutrinos decouple from the matter of the nascent neutron 
star and stream off essentially freely ( "neutrinosphere" ) is denoted by 
R^. The neutrino cooling layer and the neutrino heating region behind 
the supernova shock (at Rs) are separated by the "gain radius" Rg. 

Matter is accreted into the shock at a rate M and is partly advected 
through the gain radius into the cooling region and onto the neutron 
star. Convective overturn between the gain radius and the shock in- 
creases the efficiency of neutrino heating. Convective activity inside the 
proto-neutron star raises the neutrino luminosities and thus amplifies 
the neutrino energy deposition. 

of large-scale mixing in the stellar mantle and envelope and of fast moving Ni 
clumps. Both might indicate that macroscopic anisotropics and inhomogeneities 
were already present near the formation region of Fe group elements during the 
very early stages of the explosion. Spherically symmetric models had been sug- 
gesting for some time already that regions inside the newly formed neutron star 
and in the neutrino-heated layer around it might be convectively unstable. These 
theoretical results and the observational findings in SN 1987A were motivation 
to study stellar core collapse and supernova explosions with multi-dimensional 
simulations. 

In this article convective overturn in the neutrino-heated region around the 
collapsed stellar core is discussed concerning its effects on the neutrino energy 
deposition and its potential importance for the revival of the stalled bounce 
shock and thus for triggering the supernova explosion. Convective activity inside 
the proto-neutron star is suggested as a possibly crucial boost of the neutrino 
luminosities on a timescale of a few hundred milliseconds after core bounce. 
The first two-dimensional simulations that follow the evolution of the nascent 
neutron star for more than one second are shortly described. 
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2. Neutrino-driven explosions and convective overturn 



2.1. Neutrino heating and supernova explosions 

Figure |l| displays a sketch of the neutrino cooling and heating regions outside 
the proto-neutron star at the center. The main processes of neutrino energy 
deposition are the charged-current reactions v^ + n — > p + e~ and z^e+P — > n + e"*" 
(Bethe & Wilson 1985). The heating rate per nucleon (N) is approximately 
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where Yn and Yp are the number fractions of free neutrons and protons, respec- 
tively, ^1,^52 denotes the luminosity of or in 10^^ erg/s, the radial position 
in 10^ cm, and (e^ ^^5) the average of the squared neutrino energy measured in 
units of 15 MeV. / is the angular dilution factor of the neutrino radiation field 
(the "flux factor" , which is equal to the mean value of the cosine of the angle 
of neutrino propagation relative to the radial direction) which varies between 
about 0.25 at the neutrinosphere and 1 for radially streaming neutrinos far out. 
Using this energy deposition rate, neglecting loss due to re-emission of neutrinos, 
and assuming that the gravitational binding energy of a nucleon in the neutron 
star potential is (roughly) balanced by the sum of internal and nuclear recom- 
bination energies after accretion of the infalling matter through the shock, one 
can estimate the explosion energy to be of the order 
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AM is the heated mass, At the typical heating timescale, E'gb the (net) total 
gravitational binding energy of the overlying, outward accelerated stellar layers, 
and -Enuc the additional energy from explosive nucleosynthesis which is typically 
a few 10^*^ erg and roughly compensates Egi, for progenitors with main sequence 
masses of less than about 20 Mq. Since the gain radius, shock radius, and At 
and thus also AM depend on Ly{e1), the sensitivity of i?exp to the neutrino 
emission parameters is even stronger than suggested by Eq. 

2.2. Convective overturn in the neutrino-heated region 

Convective instabilities in the layers adjacent to the nascent neutron star are a 
natural consequence of the negative entropy gradient built up by neutrino heat- 
ing (Bethe 1990) and are seen in recent two- and three-dimensional simulations 
(Burrows et al. 1995; Herant et al. 1992, 1994; Janka & Miiller 1995, 1996; Mez- 
zacappa et al. 1997; Miller et al. 1993; Shimizu et al. 1994). Figure ^ shows 
the entropy distribution between proto-neutron star and supernova shock about 
170 ms after core bounce for one of the calculations by Janka & Miiller (1996). 
Although there is general agreement about the existence of this unstable region 
between the radius of maximum neutrino heating (which is very close outside 
the "gain radius" -Rg, i.e. the radius where neutrino cooling switches into net 
heating) and the shock position i?s, the strength of the convective overturn and 
its importance for the success of the neutrino-heating mechanism in driving the 
explosion of the star is still a matter of vivid debate. 
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Figure 2. Inhomogeneous distribution of cool and hot gas between 
the nascent neutron star (left side, middle) and the supernova shock 
front (bumpy, hemispheric discontinuity) during the first second of the 
explosion. Neutrino-heated matter rises and expands in mushroom-like 
bubbles, while cool gas flows down toward the neutron star in narrow 
streams. The radius of the semicircle is 1600 km, the shock is at about 
1400 km. 
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Figure 3. Explosion energies £^>o(i) for ID (dashed) and 2D (solid) 
simulations with different assumed v^, and Vf. luminosities (labels give 
values in lO^^erg/s) from the proto-neutron star. Below the smallest 
given luminosities the considered ISM© star does not explode in ID 
and acquires too low an expansion energy in 2D to unbind the stellar 
mantle and envelope. £'>o is defined to include the sum of internal, 
kinetic, and and gravitational energy for all zones where this sum is 
positive (the gravitational binding energies of stellar mantle and enve- 
lope and additional energy release from nuclear burning are not taken 
into account). 



The effect of convective overturn in the neutrino-heated region on the shock 
is two-fold. On the one hand, heated matter from the region close to the gain 
radius rises outward and at the same time is replaced by cool gas flowing down 
from the postshock region. Since the production reactions of neutrinos (e^ 
capture on nucleons and thermal processes) are very temperature sensitive, the 
expansion and cooling of rising plasma reduces the energy loss by reemission 
of neutrinos. Moreover, the net energy deposition by neutrinos is enhanced as 
more cool material is exposed to the large neutrino fluxes just outside the gain 
radius where the neutrino heating rate peaks (the radial dilution of the fluxes 
roughly goes as l/r^). On the other hand, hot matter floats into the postshock 
region and increases the pressure there. Thus the shock is pushed further out 
which leads to a growth of the gain region and therefore also of the net energy 
transfer from neutrinos to the stellar gas. 

2.3. Requirements for neutrino-driven explosions 

In order to get explosions by the delayed neutrino-heating mechanism, certain 
conditions need to be fulfilled. Expansion of the postshock region requires suffi- 
ciently large pressure gradients near the radius i?cut of the developing mass cut. 
If one neglects self-gravity of the gas in this region and assumes the density pro- 
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file to be a power law, p{r) oc r~" (which is well justified according to numerical 
simulations which yield a power law index of n 3; see also Bethe 1993), one 
gets P{r) oc r~^~^ for the pressure in an atmosphere near hydrostatic equilib- 
rium, and outward acceleration is maintained as long as the following condition 
for the "critical" internal energy density e holds: 

> ^ ^ - r3) 

GMp/r {n + l)[^-l) 4' 

where use was made of the relation P = (7 — The numerical value was 
obtained for 7 = 4/3 and n = 3. This condition can be converted into a criterion 
for the entropy per baryon, s. Using the thermodynamical relation for the 
entropy density normalized to the baryon density n^, s = {e+P)/{ni}T) — J2i Vi^i 
where r)i {i = n, p, e~ , e+) are the particle chemical potentials divided by the 
temperature, and assuming completely disintegrated nuclei behind the shock 
so that the number fractions of free protons and neutrons are Yp = Yg and 
Yn = I — Yg, respectively, one gets 



(i^^^,) > 15 % - In f 1.27 • 10-3 



[kB/N] . (4) 

^ut ^ ' Rc\x% 

In this approximate expression a term with a factor Yf. was dropped (its absolute 
value being usually less than 0.5 in the considered region), nucleons are assumed 
to obey Boltzmann statistics, and, normalized to representative values, M\x is 
measured in units of 1.1 Mq, in 10^ g/cm^, and r^j in 10^ cm. Inserting typical 
numbers (T k, 1.5 MeV, y„ k, 0.3, i?cut ~ 1-5 • 10^ cm), one obtains s > 15 k^/N, 
which gives an estimate of the entropy in the heating region when the star is 
going to explode. 

These requirements can be coupled to the neutrino emission of the proto- 
neutron star by the following considerations. A stalled shock is converted into 
a moving one only when the neutrino heating is strong enough to increase the 
pressure behind the shock by a sufficient amount. Considering the Rankine- 
Hugoniot relations at the shock, Bruenn (1993) derived a criterion for the heating 
rate per Tuiit mass, behind the shock that guarantees a positive postshock 
velocity (ui > 0): 

> /?3(/3-l)(7-l) ryi?, " 

Here (3 is the ratio of postshock to preshock density, (3 = p\/ Pq, 7 the adiabatic 

index of the gas (assumed to be the same in front and behind the shock), and 
r] defines the fraction of the shock radius where net heating by neutrino 
processes occurs: ry = {Rs — Rg)/Rs- uq is the preshock velocity, which is a 
fraction a (analytical and numerical calculations show that typically a ~ l/\/2) 
of the free fall velocity, uq = ayJIGMjr. Assuming a strong shock, one has 
/?= (7 + 1)7(7 — 1) which becomes /3 = 7 for 7 = 4/3. With numbers typical 
of the collapsed core of the 15 star considered by Janka & Miiller (1996), 
i?s = 200km, 77 0.4, and an interior mass M = 1.1 Mq, one finds for the 
threshold luminosities of v^. and Vf.: 

^.,52(6^,15) > 2.0^. (6) 

-^8,200 



The existence of such a threshold luminosity of the order of 2 • 10^^ erg/s 
is underlined by Fig. |^ where the explosion energy E^q as function of time 
is shown for numerical calculations of the same post-collapse model but with 
different assumed neutrino luminosities from the proto-neutron star. E^q is de- 
fined to include the sum of internal, kinetic, and gravitational energy for all 
zones where this sum is positive (the gravitational binding energies of stellar 
mantle and envelope and additional energy release from nuclear burning are not 
taken into account). For one-dimensional simulations with luminosities below 
1.9 • 10^^ erg/s we could not get explosions when the proto-neutron star was as- 
sumed static, and the threshold for the Uf, and luminosities was 2.2- 10^'^ erg/s 
when the neutron star was contracting (see Janka & Miiller 1996). The sup- 
porting effects of convective overturn between the gain radius and the shock 
described in Sect. 2.2. lead to explosions even below the threshold luminosities 
for the spherically symmetric case, to higher values of the explosion energy for 
the same neutrino luminosities, and to a faster development of the explosion. 
This can clearly be seen by comparing the solid (2D) and dashed (ID) lines in 
Fig. |. 

The results of Fig. ^ also show that the explosion energy is extremely sen- 
sitive to the neutrino luminosities and mean energies. This holds in ID as well 
as in 2D. The question can be asked why neutrino-driven explosions should be 
self-regulated. Which kind of feedback should prevent the explosion from being 
more energetic than a few times 10^^ erg? Certainly, the neutrino luminosities in 
current models can hardly power an explosion and therefore a way to overpower 
it is not easy to imagine. Nevertheless, Fig. ^ and Eq. (^) offer an answer to 
the question: When the matter in the neutrino-heated region outside the gain 
radius has absorbed roughly its gravitational binding energy from the neutrino 
fluxes, it starts to expand outward (see Eq. @) and moves away from the region 
of strongest heating. Since the onset of the explosion shuts off the re-supply of 
the heating region with cool gas, the curves in Fig. ^approach a saturation level 
as soon as the expansion gains momentum and the density in the heating region 
decreases. Thus the explosion energy depends on the strength of the neutrino 
heating, which scales with the z/g and i>e luminosities and mean energies, and it is 
limited by the amount of matter AM in the heating region and by the duration 
of the heating (see Eq. (|2[)), both of which decrease when the heating is strong 
and expansion happens fast. 

This also implies that neutrino-driven explosions can be "delayed" (up to a 
few 100 ms after core bounce) but are not "late" (after a few seconds) explosions. 
The density between the gain radius and the shock decreases with time because 
the proto-neutron star contracts and the mass infall onto the collapsed core 
declines steeply with time. Therefore the mass AM in the heated region drops 
rapidly and energetic explosions by the neutrino-heating mechanism become less 
favored at late times. 

Moreover, Fig. |3| tells us that convection is not necessary to get an explosion 
and convective overturn is no guarantee for strong explosions. Therefore one 
must suspect that neutrino-driven type-II explosions should reveal a considerable 
spread of the explosion energies, even for similar progenitor stars. Rotation in 
the stellar core, small differences of the core mass or statistical variations in 
the dynamical events that precede and accompany the explosion may lead to 
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Figure 4. Order scheme for the dependence of the post-cohapse dy- 
namics on the strength of the neutrino heating as a function of L^{e'^). 
The destiny of the star - accretion or explosion - depends on the rel- 
ative size of the timescales of neutrino heating, Tht, matter advection 
through the gain region onto the proto-neutron star, Tadvj and growth 
of convective instabihties, Tcv With a larger value of L^{e^) the heat- 
ing timescale as well decrease, the latter due to a steeper entropy 
gradient built up by the neutrino energy deposition near the gain ra- 
dius. 



some variability. The potential sensitivity to the properties and structure of the 
progenitor star, to the dynamical events during and after core collapse, and to 
the neutrino emission parameters of the proto-neutron star do not make type-II 
supernovae promising candidates for standard candles with a rather well defined 
value of the explosion energy. 

2.4. When is neutrino-driven convection crucial for an explosion? 

The role of convective overturn and its importance for the explosion can be 
further illuminated by considering the three timescales of neutrino heating, Tht, 
advection of accreted matter through the gain radius into the cooling region and 
onto the neutron star (compare Fig. Q), Tad, and growth of convective instabili- 
ties, Tcv The evolution of the shock — accretion or explosion — is determined 
by the relative sizes of these three timescales. Straightforward considerations 
show that they are of the same order and the destiny of the star is therefore a 
result of a tight competition between the different processes (see Fig. |^. 

The heating timescale is estimated from the initial entropy si, the critical 
entropy Sc (Eq. (^)), and the heating rate per nucleon (Eq. (||)) as 

^c-^i ,5^, ^c-.i fi^..(T/2MeV)/ 
^ Qt/{kBT) ^ 5kB/N {LJ2 . 1052erg/s)(e2 ' ^ > 
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With a postshock velocity of ui = uq/ (3 (7— l)yJGM / R^/ (7+I) the advection 
timescale is 

Rs-Rg _ Rg\ -^s,2oo 



where the gain radius can be determined as 
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from the requirement that the heating rate, Eq. (||), is equal to the cooling 
rate per nucleon, Q~ ~ 288(r/2MeV)^ MeV/(A^ • s), when use is made of the 
power-law behavior of the temperature according to T(r) ~ 4MeV (i?ns/^) with 
i?ns being the proto-neutron star radius (roughly equal to the neutrinosphere 
radius). The growth timescale of convective instabilities in the neutrino- heated 
region depends on the gradients of entropy and lepton number through the 
growth rate of Ledoux convection, ctl {g is the gravitational acceleration): 



dp\ d£ dYe 
dsjy.p dr \dYe)s,p dr 



'1/2 

> 50 ms 



(10) 

The numerical value is representative for those obtained in hydro dynamical sim- 
ulations (e.g., Janka & Miiller 1996). Tcv of Eq. (|lO| ) is sensitive to the detailed 
conditions between neutrinosphere (where 1^ has typically a minimum), gain 
radius (where s develops a maximum), and the shock. The neutrino heating 
timescale is shorter for larger values of the neutrino luminosity Ly and mean 
squared neutrino energy (e^), while both Tht and r^d depend strongly on the 
gain radius, Tad also on the shock position. 

In order to be a crucial help for the explosion, convective overturn in the 
neutrino-heated region must start on a sufficiently short timescale. This happens 
only in a rather narrow window of Ly{ef,) where Tcv < Tad ~ ^ht (Fig. 3). For 
smaller neutrino luminosities the heating is too weak to create a sufficiently 
large entropy maximum and the convective instability cannot develop before 
the accreted gas is advected through the gain radius (rad < t^v < Tht). In 
this case neither with nor without convective processes energetic explosions can 



occur (see also Sect. 2.3. ). For larger neutrino luminosities the neutrino heating 
is so strong (rht < Tcv ~ Tad) that expansion of the postshock layers has set in 
before the convective activity reaches a significant level. 



3. Convection inside the nascent neutron star 

Convective energy transport inside the newly formed neutron star can increase 
the neutrino luminosities considerably (Burrows 1987). This could be crucial for 
energizing the stalled supernova shock (Mayle & Wilson 1988; Wilson & Mayle 
1988, 1993; see also Sect. |]). 

Negative lepton number and entropy gradients have been seen in several 
one-dimensional (spherically symmetrical) simulations of the neutrino-cooling 
phase of nascent neutron stars (Burrows & Lattimer 1986; Burrows 1987; Keil 
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Figure 5. Profiles of the lepton fraction Yi^p = niep/ri;, (left) and of 
the entropy per nucleon, s, (right) as functions of enclosed (baryonic) 
mass for different times in a one-dimensional simulation of the neutrino 
cooling of a ~ 1.1 Mq proto-neutron star. Negative gradients of lepton 
number and entropy suggest potentially convectively unstable regions. 
Time is (roughly) measured from core bounce. 




Figure 6. Same as Fig. ^, but for a two-dimensional, hydrodynamical 
simulation which allowed to follow the development of convection. The 
plots show angularly averaged quantities in the ~ 1.1 Mq proto-neutron 
star. In regions with convective activity the gradients of l^ep and s are 
flattened. The convective layer encompasses an increasingly larger part 
of the star. 
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Figure 7. Absolute values of the velocity in the proto-neutron star 
for two instants (about 0.5s (left) and Is (right) after core bounce) as 
obtained in a two-dimensional, hydrodynamical simulation. Note that 
the neutron star has shrunk from a radius of initially about 60 km to 
little more than 20 km. The growth of the convective region can be 
seen. Typical velocities of the convective motions are several 10^ cm/s. 

h Janka 1995; Sumiyoshi et al. 1995; see also Fig. |5|) and have suggested the 
existence of regions which are potentially unstable against Ledoux convection. 
Recent two-dimensional, hydrodynamical simulations by Keil et al. (1996, 1997) 
and Keil (1997) have followed the evolution of the ~ 1.1 Mq proto-neutron star 
formed in the core collapse of a 15 Mq star for a period of more than 1.2 seconds. 
These simulations have confirmed the development of convection and its impor- 
tance for the evolution of the neutron star. They were performed with the hy- 
drodynamics code Prometheus. A general relativistic ID gravitational potential 
with Newtonian corrections for asphericities was used, $ = ^ij^ + (*^2D ~ *^id)' 
and a flux-limited (equilibrium) neutrino diffusion scheme was applied for each 
angular bin separately ("l^D"). 

3.1. Phenomenology of convection in two dimensions 

The simulations show that convectively unstable surface- near regions (i.e., around 
the neutrinosphere and below an initial density of about 10^^ g/cm^) exist only 
for a short period of a few ten milliseconds after bounce, in agreement with the 
findings of Bruenn & Mezzacappa (1994), Bruenn et al. (1995), and Mezzacappa 
et al. (1997). Due to a flat entropy profile and a negative lepton number gradi- 
ent, convection, however, also starts in a layer deeper inside the star, between 
an enclosed mass of 0.7 Mq and 0.9 Mq, at densities above several lO^^g/cm^. 
From there the convective region digs into the star and reaches the center after 
about one second (Figs. |6| and P). Convective velocities as high as 5 • 10^ cm/s 
are found (about 10-20% of the local sound speed), corresponding to kinetic en- 
ergies of up to 1-2 • lO^'' erg (Fig. ^). Because of these high velocities and rather 
flat entropy and composition profiles in the star (Fig. ^), the overshoot region is 
large (see Fig. |9|). The same is true for undershooting during the first ~ 100 ms 
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Figure 8. Convective luminosity and neutrino luminosities (dashed: 
Ljy^+Lp^, dash-dotted: L^^+Lp^+L,y^+Lp^, dotted: total) as functions 
of enclosed baryonic mass for the two-dimensional proto-neutron star 
simulation at about 500 ms after core bounce. 

after bounce. Sound waves and perturbances are generated in the layers above 
and interior to the convection zone. 

The coherence lengths of convective structures are of the order of 20-40 
degrees (in 2D!) (see Fig. |^ and coherence times are of the order of 10 ms 
which corresponds to only one or two overturns. The convective pattern is 
therefore very time-dependent and nonstationary. Convective motions lead to 
considerable variations of the composition. The lepton fraction (and thus the 
abundance of protons) shows relative fluctuations of several 10%. The entropy 
differences in rising and sinking convective bubbles are much smaller, only a few 
per cent, while temperature and density fluctuations are typically less than one 
per cent. 

The energy transport in the neutron star is dominated by neutrino diffu- 
sion near the center, whereas convective transport plays the major role in a 
thick intermediate layer where the convective activity is strongest, and radiative 
transport takes over again when the neutrino mean free path becomes large near 
the surface of the star (Fig. But even in the convective layer the convective 
energy flux is only a few times larger than the diffusive flux. This means that 
neutrino diffusion can never be neglected. 

3.2. Mechanism of convection 

There is an important consequence of this latter statement. The convective 
activity in the neutron star cannot be described and explained as Ledoux con- 
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Figure 9. Left: Convectively unstable region (corresponding to nega- 
tive values of the displayed quantity = —{g/p)Cq^ with Cql from 
Eq. ([12D) at about 500 ms after bounce according to the Quasi-Ledoux 
criterion which includes non-adiabatic and lepton-transport effects by 
neutrino diffusion. Right: Layers of Quasi-Ledoux convective instabil- 
ity (blue) and gas motions with absolute velocities larger than lO"^ cm/s 
(green) as function of time. The criterion C}^{r) = ming {Cq^{r,0)) > 
with Cql(?', ^) from Eq. ( [T^ ) was plotted. 



vection. Applying the Ledoux criterion for local instability, 

with (Tl from Eq. ( [Tol ) and replaced by the total lepton fraction Yi^p in the 
neutrino-opaque interior of the neutron star (for reasons of simplicity, Vs was 
replaced by ds/dr and Wi^p by dl^^p/dr), one finds that the convecting re- 
gion should actually be stable, despite of slightly negative entropy and lepton 
number gradients. In fact, below a critical value of the lepton fraction (e.g., 
^ep,c = 0.148 for p = 10^^ g/cm^ and T = lO.TMeV) the thermodynamical 
derivative {dp/dYiep)s,p changes sign and becomes positive because of nuclear 
and Coulomb forces in the high-density equation of state (see Bruenn & Dineva 
1996). Therefore negative lepton number gradients should stabilize against con- 
vection in this regime. However, an idealized assumption of Ledoux convection 
is not fulfilled in the situations considered here: Because of neutrino diffusion, 
energy exchange and, in particular, lepton number exchange between convec- 
tive elements and their surroundings are not negligible. Taking the neutrino 
transport effects on liep into account in a modified ^^Quasi-Ledoux criterior{\ 

V9s;(y,^^),(P) dr V^^«p/(s)(P>^ " dr J 

(12) 

(Keil 1997 and Keil et al. 1997), one determines instability exactly where the 
two-dimensional simulation reveals convective activity. In Eq. ([l^ ) the quantities 
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Figure 10. Absolute value of the gas velocity in a convecting, rotat- 
ing proto-neutron star about 750ms after bounce (left). Convection 
is suppressed near the rotation axis (vertical) and develops strongly 
only near the equatorial plane where a flat distribution of the specific 
angular momentum jz (right) has formed. 



(Ilep) and {s) mean averages over the polar angles 0, and local gradients have to 
be distinguished from gradients of angle-averaged quantities which describe the 
stellar background. The term /3iep(dllep/dr) with the empirically determined 
value Piep 1 accounts for the change of the lepton concentration along the 
path of a rising fluid element due to neutrino diffusion. Figure ^ shows that 
at about half a second after core bounce strong driving forces for convection 
occur in a narrow ring between 9 and 10 km where a steep negative gradient 
of the lepton fraction exists (see Fig. P). Further out, convective instability is 
determined only in the finger-like structures of rising, high-ljep gas. 

3.3. Accretion and rotation 

In very recent simulations, post-bounce mass accretion and rotation of the form- 
ing neutron star were included. Accretion leads to stronger convection with 
larger velocities in a more extended region. This can be explained by the steep- 
ening of lepton number and entropy gradients and the increase of the gravita- 
tional potential energy when additional matter is added onto the neutron star. 
Rotation has very interesting consequences, e.g., leads to a suppression of con- 
vective motions near the rotation axis because of a stabilizing stratification of 



the specific angular momentum (see Fig. 10), an effect which can be understood 
by applying the (first) Solberg-H0iland criterion for instabilities in rotating, self- 
gravitating bodies (Tassoul 1978): 



CsH(r,^)^4T^ + - 
x'^ ax p 



Vs 



Icp ! 



dp 



Icp 



cp 



S,P 



< . (13) 



Here, jz is the specific angular momentum of a fiuid element, which is conserved 
for axially symmetric configurations, x is the distance from the rotation axis, 
and in case of rotational equilibrium a is the sum of gravitational and centrifugal 
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Figure 11. Left: Luminosities Ljy(t) and mean energies {eu){t) of i^e 
and for a 1.1 Mq proto-neutron star without ( "ID" ; dotted) and with 
convection ("2D"; sohd). Right: Angular variations of the neutrino flux 
at different times for the 2D simulation. 



accelerations, a = SJPj p. Changes of the lepton number in rising or sinking 



convective elements due to neutrino diffusion were neglected in Eq. (13). Ledoux 
(or Quasi-Ledoux) convection can only develop where the first term is not too 
positive. In Fig. ^ fully developed convective motion is therefore constrained to 
a zone of nearly constant close to the equatorial plane. At higher latitude the 
convective velocities are much smaller, and narrow, elongated convective cells 
aligned with cylindrical regions of jz = const parallel to the rotation axis are 
indicated. 

The rotation pattern displayed in Fig. |l^ is highly differential with a ro- 
tation period of 7.3 ms at x = 22 km and of 1.6 ms at x = 0.6 km. It has 
self-consistently developed under the influence of neutrino transport and con- 
vection when the neutron star had contracted from an initial radius of about 
60 km (with a surface rotation period of 55 ms at the equator and a rotation 
period of ~ 5 ms near the center) to a final radius of approximately 22 km. Due 
to the differential nature of the rotation, the ratio of rotational kinetic energy 
to the gravitational potential energy of the star is only 0.78% in the beginning 
and a few per cent at the end after about 1 s of evolution. 

3.4. Consequences of proto-neutron star convection 

Convection inside the proto-neutron star can raise the neutrino luminosities 
within a few hundred ms after core bounce (Fig. [lT| ). In the considered collapsed 
core of a 15 Mq star Ly^ and Lp^ increase by up to 50% and the mean neutrino 
energies by about 15% at times later than 200-300 ms post bounce. This favors 
neutrino-driven explosions on timescales of a few hundred milliseconds after 
shock formation. Also, the deleptonization of the nascent neutron star is strongly 
accelerated, raising the Ve luminosities relative to the luminosities during this 
time. This helps to increase the electron fraction in the neutrino-heated 
ejecta and might solve the overproduction problem of A/" = 50 nuclei during the 
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early epochs of the explosion (Keil et al. 1996). In case of rotation, the effects 
of convection on the neutrino emission depend on the direction. Since strong 
convection occurs only close to the equatorial plane, the neutrino fluxes there 
are convectively enhanced while they are essentially unchanged near the poles. 

Anisotropic mass motions due to convection in the neutron star lead to 
gravitational wave emission and anisotropic radiation of neutrinos. The angular 
variations of the neutrino flux determined by the 2D simulations are of the 
order of 5-10% (Fig. 11). With the typical size and short coherence times of the 
convective structures, however, the global anisotropy of the neutrino emission 
from the cooling proto-neutron star is certainly less than 1% (more likely only 
0.1%, since in 3D the structures tend to be smaller) and kick velocities in excess 
of 300 km/s can definitely not be explained. 



4. Conclusions 

Neutrino-driven explosions are very sensitive to the structure of the progenitor 
star, to the details of the post-bounce dynamics, and to the properties of the 
neutrino emission from the nascent neutron star. Therefore type-II supernova 
explosions do not offer promising perspectives as standard candles. 

Convective overturn in the neutrino-heated region between gain radius and 
supernova shock can be an important help for the explosion only if the growth 
timescale of the instability is shorter than the timescale for advecting accreted 
material through the postshock region onto the proto-neutron star. This requires 
a sufficiently negative entropy gradient outside the gain radius and thus strong 
enough neutrino energy deposition. For low neutrino luminosities this is not 
the case and the star does not explode. For very large neutrino luminosities the 
heating is so fast that the postshock layers expand and an explosion develops be- 
fore convective instabilities have significantly grown. Therefore neutrino-driven 
convection is crucial for the explosion if one-dimensional models fail marginally, 
but cannot guarantee "robust" explosions independent of the size of the neutrino 
luminosities from the proto-neutron star. 

Neutrino-driven explosions are self-regulated in the sense that the explosion 
energy is limited by the amount of material in the neutrino-heated region and 
the duration of the heating. The explosion energy is of the order of (or less 
than) the binding energy of the postshock layers in the gravitational potential 
of the proto-neutron star, because the density and thus the heated mass around 
the gain radius decrease drastically as soon as the explosion gains momentum 
and the postshock gas expands away from the neutrinosphere. 

Convection inside the nascent neutron star can raise the neutrino luminosi- 
ties and the mean spectral energies of the emitted neutrinos considerably and 
thus can be a decisive boost for the neutrino heating. In addition, proto-neutron 
star convection has several advantageous properties. The first two-dimensional 
simulations show that, on the one hand, the neutrino transport out of the interior 
of the star is significantly increased only after convective activity has developed 
in a larger part of the neutron star. In the considered collapsed core of a 15 Mq 
star this requires a period of about 100-200 ms after bounce. On the other hand, 
since the deleptonization of the neutron star is accelerated, the electron neutrino 
luminosity increases relative to the antineutrino luminosity roughly between 200 
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and 400 ms after bounce. Both these findings may provide a remedy for the prob- 
lems of current supernova models, namely to develop explosions too quickly and 
thus to yield rather small neutron stars and to dramatically overproduce nuclei 
around = 50 in the ejecta (Herant et al. 1994, Burrows et al. 1995, Janka &; 
Miiller 1996). A delay of the explosion might, however, also be obtained when 
the restriction to two dimensions is dropped and the convective overturn turns 
out to be less strong in three dimensions. 

The simulations performed so far can certainly not be considered as the final 
step. Besides the need to demonstrate the existence of convective regions in the 
collapsed cores of different progenitor stars, future simulations will have to clarify 
the influence of the nuclear equation of state on the presence of convection in 
nascent neutron stars. Also, a more accurate treatment of the neutrino transport 
in combination with a state-of-the-art description of the neutrino opacities of 
the nuclear medium is needed to confirm the existence of a convective episode 
during neutron star formation and to investigate its importance for the explosion 
mechanism of type-II supernovae. Three-dimensional simulations will finally be 
required to study the influence of the direction of the turbulent cascade (from 
smaller to larger structures in 2D while it is opposite in 3D) and to settle the 
question of hot-bubble and proto-neutron star convection quantitatively. 
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